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Abstract 

Employing the concept of termination of sequence of SDS sets to describe the positive semi- 
definite property of a form, we establish a necessary and sufficient condition for deciding whether 
a given form on 7£™ is positive semi-definite or not, and show that, for a form which is (strictly) 
positive definite on TV\. , the corresponding sequence of SDS sets is positively terminating. 

The above results are exactly constructed as follows: We define the column stochastic mean 
matrix first, and then prove that if we choose countable infinite matrices from finite nxn column 
stochastic mean ones at random (repeats allowed), then the product of these infinite matrices 
will converge to a column stochastic mean matrix with rank 1. Finally, we show the proof for 
relations between termination of the sequence of SDS sets and positive semi-definite property 
of a form. 

The Maple program TSDS3, based upon these results, not only automatically prove the poly- 
nomial inequalities, but also output counter examples for those false. This method is verified to 
be very efficient and better than Polya method. 

Key words: Successive Difference Substitution; Sequence of SDS sets; Termination; Positive 
Semi-definite; Machine decision; Column stochastic mean matrix 



1. Introduction 

A traditional method which is now called difference substitution (DS) can be frequently 
found and applied to the proof of inequalities with symmetric form. However, this method 
does not involve asymmetric forms; it is also invalid for the problems failed with the first 
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substitution. To re medy that, Yang defined a general concept, difference substitution 
(lYangl .12005. 2006) , for both asymmetric and symmetric forms, and also introduced the 
method called successive difference substitution (SDS) which is the key to the studies 
later. According to these, Yang designed a heuristic Maple program SDS, which has 
been applied to prove a great many polynomial inequalities with more variables and 
higher degrees, some of which are still considered very hard with other computer algebra 
methods. 

At the 2nd Summer School on Symboli c Computation (Beijing , China, July 15-22, 
2006), Yang put forward an open problem ( Yang . July 15-22. 2 006) on the termination 
of successive difference substitution. That is, which kind of polynomials can be solved 
by SDS and which kind can not (here "can be solved" means that SDS can prove the 
polynomial inequality / > 0)? This problem, regarded as extremely difficult, has been 
spread widely. 

Furthermore, how to apply the method to deciding whether an inequality is false? 
Like traditional method, Yang's did not take it into account, which is completed in this 
paper. Thus, both in theory and application scope, the method of difference substitution 
is much stronger now than before. 

The traditional method of difference substitution chooses the upper triangular matrix 
A n as the basic substitution matrix. However, asm-> oo, A™ does not converge. Here we 
can choose the matrix K (g, > 0) as the basic substitution matrix. It is a pity that K m 
does not converge at all times when m — > oo. These make the discussion of termination 
of the sequence of SDS sets much more difficult. So , finally we select T n as the basic 
substitution matrix in order to overcome the termination problem more easily. 

A n , K and T n are respectively denoted by, 



A n 



1 1 ■•• 1 

1'-.: 

: '•■ '•■ 1 

••• 1 



K 



qi <?2 ■■■ q n 
q 2 '•• : 

••• q n 



Tn 



I 



(1) 



It is easily to prove that lim T™ exists. And T n has a better property: Let PT„ 

m—too 

be the set consisting of matrices transformed by swapping any two rows of T n , which 
comprises n\ elements. Choose countable infinite matrices from PT„ at random (repeats 
allowed), then the product of these infinite matrices will converge to a column stochastic 
matrix with rank 1. 



Main Lemma. Let fl be a set consisting of finite nxn column stochastic mean matrices 
(Dei 2.1). Choose countable infinite matrices from f2 at random (repeats allowed), then 
the product of these infinite matrices will converge to a column stochastic mean matrix 
with rank 1. 



Definition 1.1. The form (namely, homogeneous polynomials) f(xx, x 2 , ■ ■ ■ , x n ) G 7Z[xi, 
X2, • • • j x n ] 15 called positive semi-definite on 72.™ , if 

(V(x 1 ,x 2 , ■ ■ ■ ,x n ) € TV\) f(x 1 ,x 2 , ■ ■ ■ ,x n ) > 0, 



2 



where 72.™ = {(xi,X2, ■ • • ,x n )\xi > 0,x 2 > • • • ,x n > 0}. And the set of all positive 
semi-definite forms is denoted by PSD. If / > for all (x%,X2, ■ • • , x n ) ^ (0, 0, • • ■ , 0), 
then / is called a positive definite form on 72™ (briefly, a PD). 
The main results in this paper as follows: 

Theorem. 1° The form / ePSD iff the sequence of SDS sets (Def 3.4) of / is posi- 
tively terminating (Def 3.6) or not terminating. 

2° The form / £ PSD iff the sequence of SDS sets of / is negatively terminating (Def 
3.7). 

3° If the form / ePD, then its sequence of SDS sets is positively terminating. 

The paper is organized as follows. Section 2 gives the basic properties of column 
stochastic mean matrix and the proof to the main Lemma. Section 3 introduces some 
preliminary notions of difference substitution. The proof of main results and the estab- 
lishment of algorithm is presented in section 4 and the comparison with Polya theorem 
in section 5. 



2. The Column Stochastic Mean Matrix and the Proof for Main Lemma 

2.1. The column stochastic mean matrix 

A column stochastic matrix is a square matrix whose columns consist of nonnegative 
real numbers whose sum is 1 (|Horn. R. A.. Johnson. C. R. I . ll999h . The column stochastic 



matrix has a very important property, that is, let e = (1, 1, • • • , 1), and A = [au^)] is a 
column stochastic matrix, then there is eA — e. 

Lemma 2.1 Let a = (x\, x%, x n ) G 72™, where x\ + ••• + x n = I. And 

f3 = (&i, &2> • • • ) b n ) & 72™ , where the maximal element of (3 is denoted with max(/3), and 
the minimal element with min(/3) . Then there are 

1° min(/3) < X\bi + X2&2 + • • • + x n b n < max(/3). 

2° Let the subscripts that min(/3) appears on f3 are S\,S2,---, Sfe (1 < k < n) 
(that is b Sl — b S2 = ■■ ■ = b Sk — min(/3)), then min(/3) = xibi + 2:2^2 + • • ■ + x n b n iff 

Xs\ X S2 ~\~ ' ' ' ~\~ X Sk I. 

3° Let the subscripts that max(/3) appears on are ri,r 2 , • • • , r t (1 < t < n), then 

max(/3) = xibi + x 2 b 2 + • • • + x n b n iff x ri + x r2 H + x rt = 1. 

4° min(/3) = max(/3) iff all the elements of (3 are equal. 

The proof for Lemma 2.1 is obviously and is omitted. Next, there be the definition of 
column stochastic mean matrix. 

Definition 2.1 Let Q be a n x n column stochastic matrix, and /3 = (b\, 62, • • • > b n ) 
be any set of nonnegative real numbers that are not totally equal. If 

[min(/3Q), max(/3Q)] C [min(/3), max(/3)]. 

holds, then Q is a column stochastic mean matrix. 

The condition of Definition 2.1 is very strict. Consider the following two matrices, 
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Qi = 



1 

2 
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1 
2 


1 
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1 

2 





, Q2 = 


1 
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1 

2 








1 

2 


1 







1 

2 


1 

2 _ 



(2) 



in which Q\ is not a column stochastic mean matrix. Because for a vector fi = (0,0,1), 
we have [min(/3), max(/3)] = [0, 1] and [min(/3<3i), max(/3<3i)] = [0, 1], which doesn't sat- 
isfy the condition of Definition 2.1. And we can prove that Q 2 satisfies the condition of 
Definition 2.1. Thus there is a natural question, that is, what condition makes a column 
stochastic matrix become a column stochastic mean matrix? 

Lemma 2.2. The n x n column stochastic matrix A = [a^j)] is a column stochastic 
mean matrix iff for two arbitrary columns of A 

Col(A,ji) = (a (ljl) ,a (2jl) , • • • ,a (njl )) T , Co\(A,j 2 ) = (o(i,j 2 ), 0(2,j 2 ), ■ ■ ■ , a(n,j 2 )) T 
there be 

(°(l,Ji) a (l,j2)' a (2,ji) a (2,j 2 )r ■ ■ > a (n,ji)a(n,j 2 )) 7^ (0,0, ■• - ,0). 



Proof. =^:(proof by contradiction) Suppose there are two columns Col(A,j\), Col(A, j 2 ) 
of the column stochastic mean matrix A = [o(» 1-7 -)] such that 

( a (l,ji) a (l,j2)' a (2,ji) a (2,j 2 ),- ' ■ > a (n,ji) a (n,j 2 )) = (0,0, ■• - ,0). 

we might as well suppose that 

a (i,h) 7^ 0, • • • , a (t tjl) £ 0, a {t+1Jl) = 0, • • • , a {nJl) = 0, 
a (hh) = 0' ' ' ' ' °(fc,j2) = 0, 0(fc+i,j 2 ) 7^ 0, • • • , 0( n ,j 2 ) 7^ 0, 

where 1 < t < k < n. 

Let /3(f) = (V--, 1,0, ■■■,()), then [min(/3(i)), max(/3(t))] = [0,1]. By Lemma 2.1, 
t 

there be [mm((3(t)A),max(f3(t)A)} = [0,1], a contradiction. 

■<=: Let ^4 = [a^j)] be a n x n column stochastic matrix, and pick any two columns of 

it, 

Col(A,ji) = (a (1>jl) ,a (2Jl) ,---,a (nijl) ) T , Col(A,. ?2 ) = (a (1)J - 2) , a (2j2) , • • • , a ( „ j2) ) T . 
Then, by known, we have 

( a (l,ji) a (l,j2)! a (2,ji)a(2,j 2 )> ' ' ' > a (n,ji)a(n,j 2 )) 7^ (0,0, • • • ,0). 

That is, at least there is one positive integer s (1 < s < n) such that a{ s ,h) 7^ 0, d( s ,j 2 ) 7^ 0- 
Let Pi = {pi, 62, • • • j b n ) be any set of nonnegative real numbers that are not totally 
equal. And suppose that the subscripts that min(/?i) appears on /3i are Si, s 2 , • • • , Sfe (1 < 
k <n), and the subscripts that max(/?i) appears on /?i are n, r 2 , • • • , r t (1 <t < n). 
Consider the product of fii and the column ji , j 2 of A, and note 

/?2ji = kl a (lji) + ^20(2,ji) H + &n a (n,.Ji)i 
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fai 2 = &i a (ij 2 ) + b2 a (2. j2 ) + • • • + b n a(n,j 2 )- 
Suppose that faj 1 = mm(fa), faj 2 = max(fa), then by 2°, 3° of Lemma 2.1 we have 

a (si,ii) + a (s 2 ,ii) H \- a (s k , 3 i) = 1, 

a (ri,j 2 ) + a (r 2 ,j2) + • • • + a (r t .j 2 ) = 1- 

On the one hand, the matrix j4 is column stochastic, so the nonzero numbers of the 
column ji of A only appear in (ft( SlJ1 ), a (s 2 ji)''''' a (sj=ji))- Similarly, the nonzero 
numbers of the column j 2 only appear in (a( rij2 ), a(r 2 j 2 ), ' " " i a (n,j 2 ))- But there is at 
least one positive integer s (1 < s < n) such that a( s .j 1 ) 0> a (s,j 2 ) ^ 0- Thus there is at 
least one common element both in (si, S2, • • • , Sfe) and (n, r 2 , • • • , r t ). 

On the other hand, because mm(fa) ^ max(fa), then elements in (si, s 2 , • • • , Sfe, ^l, 
r 2 , • • • , r t ) are distinct. So, there is a contradiction. 

The reason for the above contradiction is that we mistakenly suppose that both 
mm(fa) and max(fa) are in faA. So mm(fa) and max(/3i) can't be in faA at one time. 
And by 1° of Lemma 2.1, thus sufficiency holds. 

For the contracting intervals, whether the length of a interval will converge to is very 
important. So, for a fixed matrix A, we should consider the maximal value of function 



mSix(XA) - min(XA) 

M A {X) = — . , , , (3) 

max(A) — min(A) 

where X = (x\, ■ ■ ■ , x n ), Xi > and max(X) — min(X) ^ 0. 

Lemma 2.3. Given a fixed column stochastic mean matrix A. The maximal value 
of function Ma(X) exists under the condition [xi > ,max(X) — min(X) 7^ 0], and it is 
less than 1. 

Proof. Firstly, suppose that min(X) = 0(otherwise use X = (xi — min(X),x 2 — 
min(X), ■ • • , x n — min(X)) instead of X). It is easily to prove that Ma(X) = Ma(X). 

Secondly, suppose that max(X) = 1 (otherwise use X = (x\/ max(X), x 2 / max(X), • • • , 
x n / max(X)) instead of X). Similarly, we can prove that Ma(X) = Ma{X). 

Thus we only consider the simplified function 



M A (X) = max(XA) - min(XA), 

where < X{ < 1 and ma,x(X) — l,min(X) = 0. 
Repeatedly use the formulas 



max(x, y) = {{x + y) + \x - y\)/2, 
mm(x, y) = {{x + y) - \x - y\)/2. 



It is easily to find that Ma(X) is continuous on the bounded closed region [0 < Xi < 
1, max(X) — l,min(X) = 0] , so the maximal value of it exists and less than 1. 

Let Ma = max(MA(X)), and we call Ma the contraction coefficient of the column 
stochastic mean matrix A. For a concrete matrix A, Ma is easy to solve. 

Example 1. Given 
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Qs = 



i 
HO 

3 2 u 

i i 1 

3 2 x 



find the contraction coefficient Mq 3 of Q3. 

Solution. Suppose that X = (x\,X2,Xs) 
the simplification of X, that is, 



there are six different possible forms after 



(0,a; 2 ,l),(l,x 2 ,0),(x 1 ,0,l),(x 1 ,l,0),(0,l,X3),(l,0,x 3 ). 



M Qs (0, x 2 , 1) = max((0, x 2 , 1)T 3 ) - min((0, x 2 , 1)T 3 ) 

= max((a; 2 + l)/3, (ar 2 + l)/2, 1) - min((a; 2 + l)/3, (x 2 + l)/2, 1) 
= 1 - (a; 2 + l)/3 < 2/3(holds when x 2 = 0). 

Analogously, it still holds for the other forms. So Mq 3 —2/3. 

Lemma 2.4. Let X = (x\, x 2 , ■ ■ ■ , x n ) is a set of nonnegative real numbers that 
are not totally equal, and A is a n x n column stochastic mean matrix. Note that the 
length of interval [a, b] is d([a, 6]), then 

d([mm(XA),max(XA)}) < M A d([mm(X),max(X)\). 

The purpose of the contraction coefficient Ma is to denote the above relation. 



2.2. Proof for the Main Lemma 

Main Lemma is not only very important for the proof of the main results but also 
very valuable itself. 

Main Lemma. Let Q, be a set that is consisted of finite n x n column stochastic mean 
matrices. Then Choose countable infinite matrices from f2 at random (repeats allowed), 
and the product of these infinite matrices will converge to a column stochastic mean 
matrix with rank 1. 

Proof. For 51 is a finite set, then there is a matrix whose contraction coefficient is 
maximal in it. Denote the maximal contraction coefficient with Mq. 
Let the column stochastic mean matrices picked in turn from f2 are 

B\ , E>2 , • • • , B m , • • • 

And let Bi — [a^j)], take the first row of Bi as the sequence of vectors fix, /3 2 , ■ ■ ■ , 

Pi = ( a (M)> a (l,2): ' ' ' 5 fl (l,n)): 
02 = (0(l,l),ffl(l,2)! ■ ■ ■ j a{l,n))B 2 , 

7 ' ' ' 7 

Pm = (0(l,l))0(l,2)j • • • , a(l,n))B2B 3 ■ ■ ■ B m , 
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If there is a positive integer k such that min(/?fc) = max(ft), then by column stochastic, 
it is clearly that the sequence {flm}m=i converges to (61, 61, • • • , b\), b\ = min(/3fc). 
Otherwise, by Definition 2.1, there is a nested closed set such that 



[min(/3i),max(/3i)] D [min(/3 2 ), max(/3 2 )] D ■■■ D [min(/3 m ), max(/3 m )] D • • • . 

And by Lemma 2.4, the length of the m th interval is less than d([min(/3i), max(/3i)])M^ 1 
Thus, when m — »■ 00, the length of the interval converges to 0. 

According to rule of intermediate sequence, the sequence of vectors /3\ , /3 2 , • • • converges 
to (61, 61, • • • , 61), where 

bi = lim min(/3 m ) = lim max(/J m ). 

Analogously, pick the prime vector /3i as the i th row of the matrix B\, then we have 
the corresponding limit vector which is denoted as 6», ■ ■ ■ , 6»). Thus 



Y[B m = 

m—l 



Then there are two corollaries about random column stochastic mean matrices. 

Corollary 2.5. For a random column stochastic mean matrix A, the sequence 
{A m }^ =1 converges to a matrix with rank 1. 

Corollary 2.6. Given a column stochastic mean matrix A. If there is a positive 
integer k such that A k is a column stochastic mean matrix, then the sequence {A m }^ =1 
converges to a matrix with rank 1. 

Corollary 2.6 is very valuable, for there are so many column stochastic matrices are 
not column stochastic mean, but there exists some positive integer k such that A k is a 
column stochastic mean matrix. The matrix Qi in (2) is just like this, and Q\ is a column 
stochastic mean matrix. 

In the next section, we will present some preliminary notions of difference substitution. 
Thus, all of the section 2 and 3 built a solid foundation for the section 4. 

3. Preliminary Notions of Difference Substitution 

n x n upper triangular matrix T n and T~ x are respectively 



h bi ■■■ bi 
b 2 b 2 ■■ ■ b 2 

bn b n • b n 



T — 



J.-1 



••• h 



1 -1 ••• 

2''. : 

1 ■■. ■■. -(n-1) 
••• n 
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Let [k\k 2 ■ • • k n ] be a permutation of 1, 2, • • • , n. P[k 1 k 2 ---k n ] = [ a (i,j)} is a n x n matrix 
for which 

a (l,fci) = 1; a (2,fe 2 ) = 1' ' ' ' i a (n,fc n ) = 1' 

and in all other positions (Permutation matrix). 

Definition 3.1. n x n square matrix B[ klk2 ... kn ] is defined as follows: 

B[ klk2 ... kn ] = P[ klk2 ... kn ]T n . 

And the set that consists of all B[ klk2 ... kn ] is denoted by PT„(there are n! elements in 
PT„). 

The standard simplex on 1Z™ is defined as follows, 

A n = {(xi,x 2 , ■ ■ ■ ,x n ) T \xi > 0, x 2 > 0, • • • ,x n > 0,xi + x 2 H h x n = 1}. 

Lemma 3.1. Given X = (xi, x 2 , • • • , x n ) T . If X £ A„, then P[ fel/s2 ... fen ]X e A„. 
Proof. The proof is very straightforward and is omitted. 

Definition 3.2. Let [kik 2 ■ ■ ■ k n ] be a permutation of 1,2, ■■■,n, and let X = 
(x\,x 2 , • • • , x n ) T . Then we set a linear mapping 

^[fcife 2 -fc„] : A" 

^ B[k ± k 2 ---k n ]X, 

Furthermore, composition of mapping ip! a ] is respectively 

V'loi] V'[a 2 ] ' ' ' V>[a m ] = VV] (^[02] ' ' ' (^[o m ] PO) ' ' ') = B [ai] B [a2] ■ ■ ■ B [am] X, 
for which '[x]' denotes a permutation of 1, 2, • • • , n. Thus, we briefly note 

^[a{\[a 2 ]-[a m ] = V^oi] VV 2 ] " " " i>[a m ] ■ (4) 

The image set of V ; [a 1 ][a 2 ]---[a m ] is denoted by Im(t/>[ 0l ][o 2 ]-[a m ])- Let Pi, • • • , P„ be points 
of A„ that arc affincly independent. The simplex with vertices Pi, • • • , P n is denoted by 
[Pi, ■■■,P n ]. Then A„ - [ei, e 2 , e„], where ej = (1, 0, • • • , 0) T , • • • , e„ = (0, 0, • • • , 1) T . 

Obviously, Im(^>[ ai ][ a2 ]...[ am ]) is the simplex [ei, e 2 , ■ ■ ■ , e n ]P[ ai ]P[ a2 ] • • • P[ 0m ]- Thus, we 
briefly note 

ImO[ ai ][a 2 ]---[a m ]) = [ei,e 2 , • • ■ , e n ]P [oi ]P [o2 ] • • • B [0m] = A n B [ai] B [a2] ■ ■ ■ B [am] . 

We use d(S) for representing the diameter of simplex S, then we have the following 
Lemma. 

Lemma 3.2. Let P\,---,P n be points of A„ that are affinely independent. The 
simplex with vertices Pi, • • • , P„ is denoted by [Pi, • • • , P n ]. Then there are 
1° 

n\ 

[Pi , Pi, ■ ■ ■ , Pn] = [J[Pl > P2, • • • 1 Pn]P[a] 7 

[a] 
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where the mark' (J ' denotes the union of the n! permutations that [a] traverses. 

[a] 

2° 

d([Pi,Fi, — ,P n ]) n-l 



d([P 1 ,P 2 ,---,P n ]B [a] ) ~ n 
Proof. (Sketch) 1° Consider the simplex [Pi, — - , P n ]-B[i, 2, •••,-«]> namely, 

rD P1+P2 Pi+P 2 + --- + P, M 
[Pl '—2 — n ] - 

Note that Pl+P2 f "' +Pi is the barycenter of simplex [Pi, P 2 , • • • , Pi], then 
1J [Pi, P2, • ■ ■ , P,i]P[ a ] is exactly the barycentic subdivision of [Pi, • • • , P n ]. 



Th e result of 2° is a elementary property of bar ycentic subdivision ( I Armstrong. M. ~ 
(jl983h : lBasu. S.. Pollack. R.. Rov. M. F. I (120061^ . 



By Lemma 3.2, we immediately have 
Lemma 3.3. 

1° 

A «= U '••UUm« 1 ] s m- b w)' 

[Om] [02] [01] 

2° 

d(A n B [ai] B [a2] ■■■B [am] ) n- 1 

d(A n ) [ n ' ■ 

Definition 3.3. When [ai], [02], • • • , [a m ] traverse all the permutations of 1, 2, • • • , n 
respectively, we define the set 

n! n! n! 

SDSM (/ ) = |J U ' ' ' U f( B M B M ■ ■ ■ B [a m ]X), 

[a m ] [o 2 ] [01] 

which is called the set of difference substitution with degree m of the form /. 

Generally, SDS^ m \f) is constructed from (nl) m polynomials. 
Now it's time to define the sequence of SDS sets. 

Definition 3.4. Given the form / £ 1Z[xi,X2, • • • , x„}, we define the sequence of sets 
{SDS (m) (/)}~ =1 as follows 

{SDS( m )(/)}^ =1 = SDS(/), SDS^(f), (5) 

We will define the termination of this sequence, which is directly related to the positive 
semi-definite property of form /. 

Let a = («i, c*2, • • • , a n ) £ N n , and let \a\ = a\ +a 2 + ■ ■ ■ +a n . Then we write a form 
/ with degree d as 

f(x u x 2 ,---,x n ) = C a x?x%''--x%>, C a £K. (6) 

\a\=d 
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Then there are some definitions on termination of {SDS*- m '(/)}^ =1 . 

Definition 3.5. The form / is called trivially positive if the coefficients C a of ev- 
ery terms x^x^ 2 • • • x" n in / are nonnegative. If /(l, 1, ■ ■ ■ , 1) < 0, then / is trivially 
negative. 

Lemma 3.5. 1° If the form / is trivially positive, then / €PSD. 
2° If the form / is trivially negative, then / g PSD. 

Definition 3.6. Given a form /, if there is a positive integer k such that every 
element of the set SDS (fe) (/) is trivially positive, the sequence of sets {SDS (m) (/)}m=i 
is called positively terminating. 

Definition 3.7. Given a form /, if there is a positive integer k and a form g such 
that g e SDS^(/) and g is trivially negative, the sequence of sets {SDS (m) (/)}^ =1 is 
called negatively terminating. 

Definition 3.8. Given a form /, the sequence of sets {SDS (m) (/)}^ =1 is neither 
positively terminating nor negatively terminating, then it is called not terminating. 

Next, let us record some properties of relations between termination of {SDS*-" 1 ' (/)}m=i 
and PSD. 

Lemma 3.6. Given the form / e TZ[xi,x 2 , • • • , x n ], and an arbitrary natural number 
m, then the following equivalence relations hold. 
1° / G PSD ^=> SDS (m) (/) C PSD. 
2° / e PSD 3g e SDS (m) (/), 5 e PSD. 
Proof. 

i° / e psd <^=> f(X) > o,vx e A„ 

^/W>0,VIe |J \J---\J(A n B [ai] B [a2 yB [am] ) 

[a m ] [a 2 ] [oi] 

f(B [ai] B [a2] ■ ■ ■ B [am] X) >0,VXe A„, VB [o<] G PT n 
-^.gePSD,V.geSDS (m) (/) 
<^ SDS (m) (.f) c PSD. 

2° getting by 1°. 



4. Proof of the Main Results 

It's time to prove main results. 

Theorem 1. Let / e lZ[x\,X2, • • • , x n ] be a form, then there are 
1° / ePSD iff the sequence of SDS sets of / is positively terminating or not terminat- 
ing. 

2° / e PSD iff the sequence of SDS sets of / is negatively terminating. 

3° If / is positive definite on 7U], then its sequence of SDS is positively terminating. 

Next we will show the proof of Theorem step by step. First prove 1° and 2°. 
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Proof of 1° and 2° in Theorem 1 

Proof. By the definitions of the termination of sequence of sets (definition 3.6, 3.7, 
3.8) we know that 1° is equivalent to 2°. So we just need to prove 2°. 

The proof of sufficiency is straightforward by the conclusion of 2° in Lemma 3.6. So 
to prove necessity, we will discuss the problem on the standard simplex A„ for / is 
homogeneous. 

Suppose that the form / 6 PSD. Then by the definition 1.1, 3X € A„ such that 
f(Xo) < 0. And we suppose that Xq is interior to A„ for such point is obviously existing. 

Since polynomials are continuous on A„, there exists a sufficient small spherical neigh- 
borhood O(X , e) of X satisfying 

VIeO(I„, £ ), f(X)<0. 

By Lemma 3.3, we can choose the simplex that include the point Xq one by one: 

A„-B[ ai ], A„B[ ai ]B[ Q2 ], A n B[ ai ]B[ a2 ] ■ ■ ■ B[ am ], ■■■ 

Then we obtain a sequence of nested closed sets 

A„-B[ Ql ] D A„B[ ai ]£?[ Q2 ] D ■■■ D A„B[ Ql ]_B[ a2 ] • ■ • B[ am ] D • • • . 

By 2° of Lemma 3.3, the diameters of these bounded closed sets are monotonically 
decreasing and tending to 0. According as the theorem of nested closed set, Xq is the 
only common point of the above sequence of closed sets. That is, when m — » +oo, there 
exists a natural number L such that 

O(X , e) D A n B [ai] B [a2] ■ ■ ■ B [ah] D ■■■ . 

Thus 

Vie A n B [ai] B [a2] ---B [aL]l f(X)<0. 
Let g = f(B [ai] B [a2] ■ ■ ■ B [i1l] X), and Y = (1, 1, • • • , 1) T . Then we have 

g(Y ) = f(B [ai] B [a2 yB [aL] Y )<0. 

That is the form f(B [ai] B [ci2] ■ ■ ■ B [ah] X) is trivially negative. Then {SDS (m) (/)}^ =1 is 
negatively terminating. 

The proof for 1° and 2° in Theorem 1 is completed. Now we'll show the proof of 3°. 

4-2. Proof of 3° in Theorem 1 

First we'll give several Lemmas for the proof of 3°. 
Lemma 4.1 Given a form of degree d 



\a\=d 



If there is g(X) = at each point X of the standard simplex A n , then g is a polynomial. 
That is, for all of the a, we have C a = 0. 

Lemma 4.1 is a well known conclusion (see Cox. P.. Little. J.. O'Shea. D. ( 19961 )). 
which can be proved by induction. 

Lemma 4.2. Given a form of degree d 
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l«|=d 

If the sequence {f m }m=i uniformly converges to on the standard simplex A„, then the 
corresponding sequence of coefficients {C( Q jTrl )}^ =1 will converge to 0. 

Proof. We only give the proof for the binary form with degree d, and the multivariate 
form can be gotten by induction. 

Let 

fm = C( di0 , m )2;f + C(o td . m )X2 + xxx 2 {C( d ^ 1 ^^ m) x d ^' 2 H V C^^i^x^ 2 ). 

Then pick (xi, X2) = (1, 0), (0, 1) and we have lim Cu m ) = 0, lim C7 dm) = 0. 
Consider the new polynomial 

/m } = /m - (C(d,0.m)xf + C\ . Am) xi) 

= x^iC^i.i^x^ 2 H h C(i id _i :m )a;2 _2 ). 

It is obviously that f$ still uniformly converges to on the standard simplex A 2 , that 
is, for all e > 0, when m is sufficient large, 

Ni^2(C( (i _i,i im )a;^ 2 H h C(i id _i jm )a;2 _2 )| < £ 

holds. 

Let ari = (1 — \f£) 1 x 2 = \fe, substituting them in the above inequality, we have 

|(1 - ^)(C (d _ lilim) (l - Vi)^ 2 + • • • + C (1 , d _ lim) (^) d - 2 )| < y/i. 

And let e — > 0, then we have C(d-i,i,m) — > 0. Analogously, let x\ = \Je, x 2 = (1 — y^e), 
then we have C(i^-i, m ) — > 0. 

Do repeated discussion, and finally we have that every sequence {(^(a^]}^- 1 con- 
verges to 0. 

Lemma 4.3. Given two forms of degree d 

fm = ^ ] £'(a,m) x l 1 x 2 2 ' ' ' x n n 1 
\a\=d 

n — r 1 -T. a lT Q 2 . . . „Q„ 

|a|=d 

If the sequence {f m }m=i uniformly converges to g on the standard simplex A„, then the 
corresponding sequence of coefficients {C( Qjm )}^ =1 will converge to C a . 

Proof of 3° in Theorem 1. We choose countable infinite matrices from PT„ at 
random: 

B [ai]> B [a 2 ], ■ ■ ■ ,B[a m ], ■■■■ 

Since PT„ is a finite set that consists of column stochastic mean matrices, by Main 
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oo 

Lemma, Yl B[a m ] * s convergent, and let 

m— 1 



b = n ^ = 

m— 1 



If we denote the degree of form / as d, X = (x\,x 2 , • • • ,x n ) T , then there always be 
f(bi, b 2 , ■ ■ ■ j 6„) > for / is positive definite on 7\L™ . Thus by the homogeneous of /, we 
know that 

00 

/( J] B [am] X) = f(BX) = f(b 1 ,b 2 , b n )( Xl +x 2 + ---+ x n ) d . (7) 

m— 1 

is trivially positive. 

Consider the infinite sequence of polynomials 

/P0> f(B [ai] X), f(B [ai] B [a2] X), /( J B [ai ]B [o2 ]---B[ am] X), (8) 

Next we will prove that there exists a positive integer mo, for to > too such that 
/( S [ai]-B[a 2 ] ' ' ' B [a m ]X) is trivially positive. 

On the one hand, for / is continuous on A„, namely, for all e > 0, there exists a 
spherical neighborhood 0{P, 5) C A n of P = (pi, 62, • • • , 6«) (If P is on the boundary of 
A„, then we take 0(P, S) n A„)) satisfying \/X u X 2 G 0(P, (5), |/(Xi) - /(X 2 )| < e. 

On the other hand, by 2° of Lemma 3.3, we know that the diameters of following 
nested closed sets are monotonically decreasing and tending to 0, 

A n D A n B [ai] D A n B [ai] B [a2] D • • • D A n B [ai] B [a2] ■ ■ ■ B [a<n] D • • • 

According as the theorem of nested closed set, these sets have only a common point. 
Obviously it is P = (61, 62, • • • , b n ). 

Putting together the above two aspects, we have that when m -4- +00, there exists a 
natural number L such that 

0(P,5)D A n B [ai] B [a2 yB [aL] D---. 

That is, when mi,m 2 > L, VX 1 e A n B [ai] B [a2] ■ ■■B [ami] (= lrn(ip [ai][a2] ... [ami] )), VX 2 € 
A n B [ai] B [a2] ■ ■ ■ B[ 0m2 ](= Im(^[oi][a 2 ]-[a m2 ])) there is 

l/pd) - /(X 2 )| < e, 

namely, VX e A„ there is 

|/(B [oi] B [oa] • • • B [ami] X) - f(B [ai] B [a2] ■ ■ ■ B [am2] X)\ < e 

So the sequence of polynomials (8) uniformly converges to f(BX) on A n . By Lemma 
4.3, the sequence of coefficients of the same monomials in the sequence (8) will con- 
verge to the coefficients of the same monomial of f(BX). And the coefficients of f(BX) 



bi bi ■■■ bi 
b 2 b 2 ■■ ■ b 2 

bn b n • b n 
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are all positive (see (7)), so there exists a positive integer m , for m > m such that 
f{B[ ai ]B[ a2 \ ■ ■ ■ B[a m] X) is trivially positive. 

For the infinite sequence (8) is arbitrary, the sequence of sets {SDS^ m ^(/)}^ =1 is pos- 
itively terminating. 

4-3. Corollary, algorithm and other consideration 

It's well-known to us that if the binary form / is a positive semi-definite form with no 
square factors, then / is definitely a positive definite form. 

Corollary 4.4 The binary form / has no square factors, then / is positive semi- 
definite iff the sequence of sets of / is positively terminating. And / is not a semi-definite 
form iff the sequence of sets of / is negatively terminating. 

Thus we obtain the following algorithm by the theorem 1, which is used to decide the 
positive semi-definite form. 

Algorithm(NEWTSDS) 

Input: form / <G Q[xi,x 2 , ■ ■ ■ ,x n ]. 
Output: the form / G PSD, or ,/ePSD. 
Tl: Let F = {/}. 

T2: For each element in F Compute [j SDS(g). 

g<BF 

Let Temp = (J SDS(g) — {trivially positive polynomials in \J SDS(g)}. 

g€F gEF 

T21: If Temp is null, then output / ePSD and terminate. 

T22: If there are trivially negative polynomials in Temp, then output 

/ G PSD and terminate. 
T23: Else, let F = Temp, go to step T2. 

According to the above algorithm, we designed a Maple program called TSDS3. It 
mainly has two orders, TSDS and NEWTSDS. TSDS chooses A n as basic substitution 
matrix, and NEWTSDS chooses T n as basic substitution matrix. 

To the order TSDS, there are some PD forms making the program do not terminate. 
For example, 

x 2 

f = (3x 1 +x 2 -x 3 ) 2 + ^-. 

And to the order NEWTSDS, by Theorem 1, we confirm that the program will termi- 
nate if the input form / itself is not PSD, and judge that / > does not hold. Otherwise, 
if the input form / itself is PSD, then the program either terminate(that is, judge that 
/ > holds) or does not(namely, can't decide whether / > holds). However, as the 
3° in Theorem 1 says, the algorithm NEWTSDS is complete for the positive definite 
form. And if we enter the subprogram of eliminating square factors, then the algorithm 
NEWTSDS is also complete for the decision of binary form. 

Then it follows a question naturally: What is the necessary and sufficient condition to 
determine that a {SDS^ m ^(/)}^ =1 is not terminating for a general form /? 

We can generalize the Theorem 1 as follows 
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Let R\, i?2) • • • , Rk be n x n column stochastic mean matrices, such that 

k 

[Pl,P2, ■ ■ ■ , Pn] = |J [Pi, P2, ■ ■ ■ , Pn]fli, 

where [Pi, P2, ■ • • , -P n ] is arbitrarily a simplex in A„. And the set SDS^ m ^(/) is definited 
as 

k k 

SDS (m \f)= \J ... (J /:/,', •••/.', AK 

2i=l i m -l 

Then the results of Theorem 1 still hold. 

In addition, from the aspect of the method of partial cylindrical algebraic decomposi- 
tion (PCAD) (jCollins, G. E.. Hong. H~l . Il99lh . it is easily to prove that: 

Lemma 4.5. Given the form / £ 1Z[xi,X2, • • • , x n ]. There exists a positive integer m 
such that the point set 



W [ Q 2] [ Q l] 



has at least one sample point in each cell (here all the lower dimensional cells are dis- 
carded) . 

Lemma 4.5 gives another way for the completion of the algorithm NEWTSDS. That 
is, to compute the upper bound of to, which is a fine subject for further study. 



5. Comparison between Theorem 1 and Polya's Theorem and Applying ex- 
amples 

Theorem (Polya's) If a form / £ lZ[xi,X2, ■ ■ ■ ,x n ] is positive on 72.™ when x\ + 
X2 H + x n > 0, then for sufficiently large N all the coefficients of 

(xi + x 2 H h x n ) N f 

are positive. 

Polya's Theorem appeared in 1928 (in German) a nd is also in ^Inequalities^ by 
Hardy. G. H.. Littlewood. J. E.. and Polva. G. I ( 1952). From th en on, this elegant and 
beautiful result has many applications. In 1946, Habicht ( 1940| ) skillfully used Polya's 
Theorem to gi ve a constructive proof of a pos itive definite case of Hilbert's 17th Problem. 
And recently JSchweighofer. M. I (|2002l . \2004 used Polya's Theorem to give an algorithm 
proof of Schmitdgen's Positivstellensatz. Furthermore, the quantitative form of Polya's 
Theo rem has been used in the study of some other problems ( de Klerk. E.. Pasechnik. D. I . 
20021 ). and the study of quantitative problems themselves obtains very good results 



( Powers. V.. Reznick. B. I . 12001 . 20061) . We hope sincerely that the results in this pa- 
per will do some similarly contributions on the above fields. 

Next we will show the comparison between Theorem 1 and Polya's Theorem in several 
following aspects: 

Table 1 
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The given classes 


Polya's Theorem 


Theorem 1 


PD 


Completely decide 


Completely decide 


PSD with boundary zero 


Partially apply 


Partially apply 


PSD with interior zero 


Can't decide 


Partially apply 


not PSD 


Can't decide 


Completely decide 



From the above table, we know that Polya's Theorem holds at item 1, but fail at item 
3 and 4, while Theorem 1 in this paper is complete at item 1 and 4, and partial applicable 
at item 2 and 3. Next we will give some practical examples to show this difference. 

Example 2.(PD, A-G) Let (x\, X2, ■ ■ ■ , £3) € 7?.", and e > is small enough. 
Consider the positive definite form 



i=i 2=1 

fn the book o f lHardv. G. H.. Littlewood. J. E.. and Polva. G. I (Il952lh this example is 
analyzed by Polya's method, where when TV equals " 2e > we nave a trivially positive 
polynomial. 

Analogously, the Hurwitz identity relation is recorded in the same book, 



i=l i=l 



x 2 2 )( x l - X2)X 3 



x 2 3 )( X 1 - x 2) x Z x i H ), 



where the notion J^! denotes the sum of all permutations of the variables. From the 
Hurwitz identity relation we easily get that the elements of the set SDS(/i) is trivially 
positive(the basic substitution matrix can be any one of A n and T n ). 



Example 3.(PSD, Schur) (lHardv. G. H.. Littlewood. J. E.. and Polva. G. I (Il952h . 
Mitrinovic. D. S.. Vaic. P. M. I (Il970ll 'l Let x, y, z are nonnegative numbers and t is a 
real number, and let x 7^ of x° as usual. Then 



.9 



x t (x - 



y)(x - z) + y*(y - x)(y - z) + z\z - x)(z - y) > 



(9) 



with equality if and only if x = y = z or all permutations of x = 0, y = z. 

For g has interior zero (1,1,1), then it is not suitable for Polya's Theorem. And 
we'll give the analysis of Theorem 1. Using the basic substitution matrix A%( or T3), we 
compute SDS(g). By symmetry, we only compute g(A-^X), X — (x,y,z) T . 

g(A 3 X) = (x + y + zfx 2 + ((x + y + zf - (y + zf + z*)yx + z l y 2 . 



1G 



Note that (x + y + z) 1 — (y + z) 1 + z t > obviously holds, so inequality (9) holds, with 
equality if and only if 

((x + y + zf - (y + zf + z v )yx = 
zV = 

(x + y + zfx 2 = 

Solve that y = 0, x — 0, and z is arbitrary, or x = 0, z = 0, and y is arbitrary. And then 
return to original variables by multiplying A% , thus we get the condition of equality. 
That is, x — y — z or all permutations of z = 0, y = x. 



Example 4 (not PSD, Vase conjecture ) (|Yang. L.. Xia. B. C. 1 . 120081) Try to decide 
whether the inequality holds. 



<2l — &2 &2 — 03 03 — (24 CJ4 — 05 05 — (3,6 Cl6 — ai 



> 0, 



a 2 + 03 a 3 + 04 04 + a 5 05 + a 6 a 6 + ai ai + a 2 



in which a.; > 0(i = 1, • • • , 6). 

Take off denominators of the left polynomial, and then denote the new polynomial by 
/. Utilizing the program TSDS and executing order tsds(f), we have a counterexample: 
ai = 84, a 2 = 7,03 = 79,a4 — 5,05 = 76, <26 = 1. Analogously, it is easy to prove that 
the inequalities like above in 3 ,4,5,7 variab l es (in case of 7 needs longer time) still holds 



There are more examples in ([Yang 1,12005, 2006l |julv 15-22. 2006l:lYang. L.. Xia. B. C 



120081 ) and on the web |http://www.irgoc.org/bbs| /dispbbs.asp? boardid=12 &id=3105 
&page= &star=10, we do not list them one by one. 
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